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Abstract 

We discuss a cosmological solution of the system which was originally 
introduced by Calogero and is today popularly known as "Nordstrom- 
Vlasov system". 

Although the model is unphysical, its cosmological solution results 
interesting for the same reasons for which the Nordstrom- Vlasov system 
was originally introduced in the framework of galactic dynamics. In fact, 
it represents a theoretical laboratory where one can rigorously study some 
problems, like the importance of the gravitational waves in the dynamics, 
which at the present time are not well understood within the physical 
model of the Einstein- Vlasov system. 

As the cosmology of the Nordstrom- Vlasov system is founded on a 
scalar field, a better understanding of the system is important also in 
the framework of the Dark Energy problem. In fact, various attempts to 
achieve Dark Energy by using scalar fields are present in the literature. 

In the solution an analytical expression for the time dependence of 
the cosmological evolution of the Nordstrom's scalar field is also released. 
This analytical expression is unknown in the literature of the Nordstrom- 
Vlasov system. 

Based on their importance, the propagation of gravitational waves in 
the Nordstrom- Vlasov system and their effects on test masses are also 
discussed. 

Keywords: Nordstrom- Vlasov system; gravitational waves, 
Dark Energy. 
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1 Introduction 



Historically, the relativistic scalar theory of gravitation introduced in 1912-13 
by the Finnish physicist Gunnar Nordstrom [T] has been the first metric theory 
of gravity. In fact, it was derived three years before Einstein's general relativity 
[2J. Although Nordstrom's scalar theory of gravity is today ruled out by obser- 
vations, it results an interesting example of non-minimal coupling of matter to 
gravity and shares with general relativity the unique property to embody the 
strong equivalence principle [3]- It is also considered an useful toy model to test 
analytical [4] or numerical [5] analyses on the two-body problem. 

The interesting idea to couple Nordstrom's scalar theory of gravity with 
the Vlasov's kinetic model for a collisionless ensemble of particles interacting 
through the gravitational force is due to Calogero, who analyzed the "Nordstrom- 
Vlasov system" in the framework of galactic dynamics [6] . Calogero and collab- 
orators further discussed the Nordstrom- Vlasov system and its solutions in refs. 

mm- 

In this paper we discuss a cosmological solution of the Nordstrom- Vlasov sys- 
tem. The cosmological framework is interesting for the same reasons for which 
the system was originally introduced in the framework of galactic dynamics [B]. 
It is considered an interesting theoretical laboratory where one can rigorously 
study problems which at the present time are not well understood within the 
physical model of the Einstein- Vlasov system [6J . An example is the importance 
of the gravitational waves in the dynamics [5J . 

We recall that the cosmology of the Nordstrom- Vlasov system is founded on 
a scalar field. Hence, a better understanding of the system is important also 
in the framework of the Dark Energy problem as various attempts to achieve 
Dark Energy by using scalar fields are present in the literature see a d 
references within. 

The solution shows also the cosmological evolution of Nordstrom's scalar 
field which is unknown in the literature of the Nordstrom- Vlasov system. 

Based on their importance, the propagation of gravitational waves in the 
Nordstrom- Vlasov system and their effects on test masses are also discussed in 
the last Section of this paper. 

2 The Nordstrom- Vlasov system 

Let us clarify an important issue. One has to distinguish between the standard 
kinetic theory, which is based on the Boltzmann equation, and the Vlasov's 
framework and its applications, i.e. the electromagnetic Vlasov-Poisson sys- 
tem, the Vlasov-Poisson non-relativistic system and the relativistic Nordstrom- 
Vlasov and Einstein- Vlasov systems. In a kinetic description the time evolution 
of the system is determined by the interactions between the particles, which 
depend on the physical situation [131 Q3] . For instance, the driving mechanism 
for the time evolution of a neutral gas is the collision between particles (the 
Boltzmann equation) [T31 Q3] . For a plasma the interaction is through the elec- 
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tromagnetic field produced by the charges (the Vlasov-Maxwell system), and in 
astrophysics the interaction between collisionless particles is gravitational (the 
Vlasov-Poisson system and the Nordstrom- Vlasov and Einstein- Vlasov systems) 
[T31 [H]. Thus, an important difference between the standard kinetic theory and 
the Vlasov's framework is that in the first case there are collisions between par- 
ticles while in the Vlasov's framework particles interact electromagnetically or 
gravitationally only and they are collisionless |1 31 114| . 

If one wants to satisfy the condition demanding that the particles make up 
an ensemble with no collisions in the spacetime, the particle density must be a 
solution of the Vlasov equation |14l [T5] 

*/ + ^W-rj^<w = o. (i) 

In this paper we work with 8ttG = 1, c = 1 and h = 1 , while the sign conventions 
for the line element, which generate the sign conventions for the Riemann/Ricci 
tensors, are (— , +, +, +). Latin indices run from 1 to 3, and Greek from to 3 
and x° = t. Here represent the usual connections, / is the particle density 
and p° is determined by p a (a = 1, 2, 3) according to the relation [TH [j~5] 

9^P V = -1, (2) 

which expresses the condition that the four momentum p M lies on the mass shell 
of the metric. 

We recall that, in general, the Vlasov-Poisson system is [THUS] 
dtf + v- Vx/ - VxU ■ v J = 



AU = 4np (3) 

p(t,x) = J dvf(t,x,v), 

where t denotes the time and x and v the position and the velocity of the 
particles. The function U = U(t, x) is the average Newtonian potential gener- 
ated by the particles. This system represents the non-relativistic kinetic model 
for an ensemble of particles with no collisions, which interacts through the grav- 
itational forces that they generate collectively [TH I15|. Thus, one can use such 
a system to describe the motion of galaxies within the Universe, thought of as 
pointlike particles, when the relativistic effects are negligible [T?]. In this ap- 
proach, the function f(t,x,v) in the Vlasov-Poisson system © is non-negative 
and gives the density on phase space of the galaxies within the Universe [14] . If 
the relativistic effects become not negligible, the motion of the particles has to 
be described by the Einstein- Vlasov system [6]. In this relativistic system, the 
dynamics of the matter is still described by a non- negative scalar function /, 
but the unknown of the field equations is now the metric of the space-time [B]. 

The Einstein- Vlasov system is much more complicated than the Vlasov-Poisson 
system [6j. The greatest difficulties are due to the hyperbolic and highly non- 
linear character of the Einstein field equations, to the equivalence of all the 
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coordinate systems in general relativity and to the issue that in this theory of 
gravitation the space-time is itself part of the solution of the equations [6] . In 
fact, obtaining global solutions of the Einstein field equations coupled to any 
kind of matter is extremely difficult |16l IT?] . 

For these reasons Calogero [6] proposed a simplified model, the Nordstrom- 
Vlasov system, in which the dynamics of the matter is still described by the 
Vlasov equation but where the gravitational forces between the particles are 
now supposed to be mediated by Nordstrom's scalar theory of gravity PQ. In 
such a theory the gravitational forces are mediated by a scalar field <j> an d the 
effect of such forces is to induce a curvature in the space-time [UG>]. Moreover, 
the scalar field modifies the otherwise flat metric only by a rescaling [TJ [S]. 
Hence, the spacetime of the model is given by the conformally flat metric [TJ [5] 

ds 2 = A 2 ((j))(dt 2 -dz 2 -dx 2 -dy 2 ) . (4) 
If one requests the property of scale invariance the function A results set as 

A{4>) = exp(^) (5) 
(except a constant which is fixed equal to one), see [5] for details. 



3 Cosmological solution 

By assuming a time dependence of the Nordstrom's scalar field <f> = <fi(t). the line 
element ^ results analogous to the cosmological Friedmann-Robertson- Walker 
(FRW) line element of the standard homogeneous, isotropic and flat Universe 

muz] 

ds 2 = A 2 [<j>(t)] (dt 2 - dz 2 - dx 2 - dy 2 ) , (6) 
and the function A [<f>(t)] results to be the scale factor. 

By using the classical transformation from conformal time to synchronous 

time mrjun] 

dt 

the line element in spherical coordinates, becomes 

ds 2 = dt 2 - A 2 [<j>(t)\ (dr 2 + r 2 (d0 2 + sin 2 9dip 2 )). (8) 
The metric tensor has the form |141 [T6l [T7] 

1 

9^ =1 I , (9) 



In 



where j mn = A 2 [<f>(t)] . 
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Let us define [33] Xmn = ^7 mn . The Einstein field equations in the syn- 
chronous frame are |141 [TBI [T7] 

^ = "5^x2 - ^xU 1 = (t ° - ir) (io) 

K = l(xl, b -xl, a )=T2 (11) 

< = -i^ - ^§- t (VlX b a ) = Tl ^T, (12) 

where P% is the Ricci tensor in 3 dimensions |14[ [TBI IT7] . 
On the other hand, the Einstein field equations in the Vlasov's framework 
are [HOI] 

Gfiu = y f(t,x,p)p tl p l/ 8(p 2 + m 2 )d 4 p, (13) 

where m is the mass of a particle (galaxy). 

Following [14] we can split the function f(t,x,p) into a couple of equations 
for x,p) and /_(£, £,p) which are constructed by reducing f(t, x,p) respec- 
tively on the "upper" half and on the "lower" half of the mass shell. Eq. JTJ 
becomes [14] 

If d 1 d~i nr d \ 

dtf± = - — 7">« ^ - 5 S—PnVr ^~ /±- (14) 

Eq. (|14p can be interpreted in Hamiltonian terms |14| 

P° ± d t f ± = {HJ ± } : (15) 
where the Hamiltonian function is |14| 

H EE ^ mn Pm Pn- (16) 

One can calculate the components of energy-momentum tensor T M „ in the 
approximation which considers galaxies like massless particles P3] (m = in 
Eq. Q3J) 



T 1 / f+ + f- . hl+Pl , gj ,3„ M7 s 

00 A 3 [0(i)]r 2 sin(9,/ A [0(t)] V r 2 r 2 sin<T^ 1 ; 

1 f a,^ (/+ + /-) „ „ * 



Tm " = TFfTTTvi 2 ■ a / r^— , ' , PmPnd p (18) 

A 3 [0(i)] sin 6> J I pi+pi p \ 

„2 I „2 „;„ a 
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The Einstein field equations (fTU|). (JTTJ) and (jTSJ give two independent dy- 
namic equations which can be written down in terms of the scale factor A(<p): 



1 [m - -am - 2 + I (/+(s) + f - {s))d ^ (21) 

where [T3] 

^Pl + 4 + ^T- n - ( 22 ) 
r 2 sin 9 

By introducing the dimensionless variables r and t |14| we put 

A[<f>(t)] = A(<h)r 
t = A(<fo)t 

(23) 

dr 



j = iA 2 (0 o )p o , 

where A(<po) = A [(f)(to)] is the present value of the scale factor of the Universe 
being </>o = cj>(to) the present value of the Nordstrom's scalar field. 
Eq. O becomes [T4] 

f 2 = -l + i 



(24) 

The solution of the system ([24]) is [T4] 



Eo = l- 



£(*) = Vi"(<- VT^l) 2 (25) 

if j > 1. Notice that expression (|25[) becomes imaginary for certain values 
of the parameter t. Such a parameter represents the cosmic time normalized by 
the present value of the scale factor A(<fio), see eq. (|23|) . By choosing the origin 
of the cosmic time at the present era of the cosmological evolution, i.e. to = 0, 
expression (|25p is always real if j > 1 . 

Returning to the (t, A [<fi(i)]) variables we get: 



A[m = m ^^- ( ^^^- ir . (26) 

The today's observed Hubble radius and the today's observed density of the 
Universe are respectively [18] A(<j>o) > 10 28 cm and po ps 10" 
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If one inserts these data in Eq. (|2"6")l a singularity at a time 
t s ps — 10 10 seconds and a value for the today's theoretical Hubble constant 
H = 4|^y « lO-^cm" 1 are obtained. 

Although the discussed model is unphysical, we note that, even under the 
assumption to neglect the baryon mass of the galaxies, the results look reason- 
able as they are of the same order of magnitude of the standard cosmological 
model QJ]. 

By using eq. (JSJ) , eq. (|2"6")l permits to obtain the cosmological evolution of 
the Nordstrom's scalar field 



<j>(t) = O + i In 



po cxp(20 o ) / t / p exp(2</> ) 
3 I cxp(^o) 



(27) 



We emphasize that this is the first time that an analytical expression for 
the time dependence of the Nordstrom's scalar field is released. In fact, the 
expression (|27|) is unknown in the literature of the Nordstrom- Vlasov system. 



4 Propagation of gravitational waves and motion 
of test masses in the Nordstrom- Vlasov system 

Let us linearize the line element (j4]) . Putting <f> = 4>* + ^ with h -C 1 and using 
([5]), eq. becomes 

ds 2 = exp (20, + h) (dt 2 - dz 2 - dx 2 - dy 2 ) = 

(28) 

= exp(2</>») exp(h) [dt 2 — dz 2 — dx 2 — dy 2 ) . 

We can also normalize the line element (|28|) to cxp(2</>*) without loss of 
generality. Hence, to first order in h, 

ds 2 = (l + h) (dt 2 ~ dz 2 ~ dx 2 - dy 2 ) . (29) 

As we are interested to the propagation of gravitational waves, we consider 
the linearized theory on vacuum by setting equal to zero the right hand side 
of the Einstein field equation (fT3|) . Following [16], to first order in h, we call 
Rfivpa , Rfiv and R the linearized quantity which correspond to R^ piJ , R^u and 
R. By inserting the line element (|29|) in the linearized field equations 



Rnu - \i\?u = 0, (30) 
we obtain only one independent equation, the d'alembertian 

□/i = 0. (31) 
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The solution of eq. ([51]) is the plane wave 

h = a{lc) exp(ik a x a ) + c.c. (32) 

which represents a scalar massless mode. Considering a wave propagating 
in the positive z direction, the line element (|29[) reads 

ds 2 = [1 + h(t - z)} (dt 2 - dz 2 - dx 2 - dy 2 ) . (33) 

Now, let us discuss the effect of these gravitational waves on two test masses 
mi and m-i. We will use the coordinate system in which the space-time is locally 
flat and the distance between any two points is given simply by the difference 
in their coordinates in the sense of Newtonian physics |16| . This frame is the 
proper reference frame of a local observer, located for example in the position O 
of the first test mass m\ . In this frame gravitational waves manifest them-self 
by exerting tidal forces on the two test masses [IB]. A detailed analysis of the 
frame of the local observer is given in Sect. 13.6 of [TB]. Here, only the more 
important features of this frame are recalled JTB] : 

the time coordinate xq is the proper time of the observer O; 

spatial axes are centered in O; 

in the special case of zero acceleration and zero rotation the spatial coor- 
dinates Xj Eire the proper distances along the axes and the frame of the local 
observer reduces to a local Lorentz frame. In this case the line element reads 

m 

ds 2 = -(dx ) 2 + .),,<!. r'd.r* + 0(\x j \ 2 )dx a dx 13 . (34) 

The effect of the gravitational wave on test masses is described by the equa- 
tion ng 

x l = —R l oko x k , (35) 

which is the equation for geodesic deviation in this frame. 

Thus, to study the effect of the gravitational wave on test masses, R l oko has 
to be computed in the proper reference frame of the local observer. But, because 
the linearized Riemann tensor R^ vpa is invariant under gauge transformations 
|16| . it can be directly computed from eq. (|33p. 

From [16] one gets 

Rp.v P a = ^{d[j,dph alJ + d u d a h^p - dadph^ - d^h a fj}, (36) 
that, in the case eq. ([33]), as a sole polarization is present, begins 

^o 7 o = \{d a d h^ + dod^hSg - d a d 1 h Voo - d d h5«}; (37) 

the different elements are (only the non zero ones will be written down 
explicitly): 



8 



d a d hr, { 



0-) 



dod-yhSy = 



dfh for a = 7 = 

—d z dth for a = 3;7 = 
d 2 h for a = 7 = 

dtd z h for a = 0; 7 = 3 



(38) 
(39) 



d a dyhrioo = d a d 1 h = < 



-dfh for a = 7 = 
d 2 h for a = 7 = 3 
dtd z h for a = 0; 7 = 3 
<9 z 9 t /i /or a = 3; 7 = 
dodohS" = —d 2 h for a = 7 . 



(40) 



(41) 



Now, putting these results in eq. (|37|) one gets 



pi _ _ 1 'i 
-"•010 — 2 



R2 
rt 010 



■±/l 



t? 3 — in/) 

-"•030 ~~ 2 



The third of eqs. (05]) together with eq. ([31]) gives R 



030 



(42) 



0. Hence, we 



have shown that gravitational waves in the Nordstrom's theory of gravity are 
transverse to the direction of propagation. In fact, using eq. (|35[) one gets 



-hx, 



\hy. 



(43) 



(44) 



Equations (f4"3"|) and (f4"4"|) give the tidal acceleration of the test mass caused 
by the gravitational wave respectly in the x direction and in the y direction. 
Equivalently, one can say that there is a gravitational potential |16| 



V(l>,t) = --h(t)[x 2 +y 2 ], 



(45) 



which generates the tidal forces, and that the motion of the test mass is 
governed by the Newtonian equation |16) 
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(46) 
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The solution of eqs. (|4"3")) and (|4"4")) can be found by using the perturbation 
method [16]. To first order in h, the displacement of the test mass due by the 
gravitational wave are given by 

Sx(t) = l -x h{t) (47) 

and 

where xq and j/o are the initial (unperturbed) 

m 2 . 

Conclusion remarks 

In this paper a cosmological solution of the Nordstrom- Vlasov system has been 
discussed by releasing an analytical expression for the time dependence of the 
Nordstrom's scalar field which is unknown in the literature. Although the model 
has to be considered unphysical, this cosmological solution results interesting 
for the same reasons for which the Nordstrom- Vlasov system was originally 
introduced in the framework of galactic dynamics in [6]. In fact, it represents a 
theoretical laboratory where one can rigorously study some problems, like the 
importance of the gravitational waves in the dynamics, which at the present 
time are not well understood within the physical model of the Einstein- Vlasov 
system. 

Based on their importance, the propagation of gravitational waves in the 
Nordstrom- Vlasov system and their effects on test masses have also been dis- 
cussed in the last Section of this paper. 

The cosmology of the Nordstrom- Vlasov system is founded on a scalar field. 
Thus, a better understanding of the system is important also in the framework 
of the Dark Energy problem. On this issue, we recall that various attempts to 
achieve Dark Energy by using scalar fields are present in the literature |19]-[23j 



(48) 

coordinates of the test mass 
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